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Abstract 



This article outlines the search for an exact general relativistic description of the exterior 
(vacuum) gravitational field of a rotating spheroidal black hole surrounded by a realistic 
axially symmetric disc of matter. The problem of multi-body stationary spacetimes is first 
J> ■ exposed from the perspective of the relativity theory (section ||) and astrophysics (section ||) , 

listing the basic methods employed and results obtained. Then (in section 3) basic formulas 
for stationary axisymmetric solutions are summarized. Sections ^ and ^ review what we 
have learnt with Miroslav Zacek and Tomas Zellerin about certain static and stationary 
■ situations recently. Concluding remarks are given in section ^. Although the survey part 

£S| \ is quite general, the list of references cannot be complete. Our main desideratum was the 

informative value rather than originality — novelties have been preferred, mainly reviews 
' and those with detailed introductions. 

cr 

1 The fields of multi-body systems involving black holes 

The subject of self-gravitating sources around rotating black holes is interesting in several re- 
spects, relevant from the point of view of the relativity theory itself as well as in the astrophysical 
context. 

First, due to the non-linearity of Einstein's equations, the field of a multi-body system is a 
traditional challenge where one typically does not manage with a simple superposition.]] On the 
first post-Newtonian level, the "celestial mechanics" of gravitationally interacting bodies can be 
kept linear |89[ , |90[| , but in the strong-field region the interaction may bring surprising features 
that have only been described in a very few cases yet (for a two-body problem — today at the 
centre of attention because of the expected gravitational waves from colliding compact binaries, 
"the current state of art is the third post-Newtonian approximation" |91[|). 

The second point in which the problem embodies the essence of general relativity is the 
effect of inertial frame dragging due to the rotation of the sources. Contrary to the Newtonian 
treatment which does not discriminate between static and stationary situation, the field is now 
determined not only by mass-energy configuration, but also by its motion within the bodies. 
The inertial space can be imagined as a viscous fluid mixed by the sources. In today's "gravito- 
electromagnetic" language, gravity has not only an electric component, generated by the mass, 



but also a magnetic one, generated by mass currents; see [245, 103, 215, |298|, 147, 79, 202, 201 



Even the problem of an "isolated" object is very difficult (mainly if its motion must be found as well, not 
mentioning the interior field) unless one makes some simple assu mptions about its multipole structure; see e.g. 
| p8| , |99| |l7| for a treatment of bodies with matter interior, [297] for that also valid for singular bodies such as 
black holes, and |40| for the case of point-like particles. 
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|163| ] and references therein. (The long-time effort to measure the gravitomagnetic effect of the 
Earth is just culminating ff8| ; cf. [f238| , p39(| .) 

The third point is that our setting involves a black hole, perhaps the most bizarre of the 
new predictions of general relativity]^] Within and around it, the deviations from Newtonian 
theory become dominant, in particular the above mentioned implications of non-linearity and 
rotation reveal themselves prominently. Perhaps the most extreme example (of an exact space- 
time) involving all the three aspects^] is a "double Kerr(-Newman)" solution for an aligned pair 
of rotating ultracompact objects [96, |210| , 37, |211| ] (cf. [ j!76| ] where an approximate method for 



axially symmetric superpositions of rotating black holes was proposed as a starting point of 
computations of black hole collisions). 

Relativistic spacetimes are found in three ways: by numerical solution of Einstein equa- 
tions, by perturbation of previously known spacetimes, and by exact analytical solution of 
Einstein equations. Let us mention what the above approaches told us about stationary ax- 
isymmetric spacetimes describing rotating black holes with additional matter (ring, disc or 
torus) .0 

1.1 Numerical solution 

Although computers have been naturally employed in relativity for decades, in algebra as well 
as in numerics, it is only quite recently that they embarked on numerical solution of Einstein's 
equations in the most interesting cases involving black holes. Teukolsky j]295| l lists today's peak 
parallel-type hardware and software and the appropriate, hyperbolic formulation of the field 
equations as "three reasons why we are on the verge of important advances in the computer 
solution of Einstein's equations". Whereas analytical prospects are restricted, the present-day 
computational facilities can handle almost any situation, at least for certain period of time. On 
the other hand, it is never sure whether a given numerical solution represents a typical or just 
a marginal case. Due to this lack of generality, it is difficult to discern, analyse and interpret 
different classes of solutions within boundless ranges of possibilities. However, numerical solu- 
tions can provide explicit examples of spacetimes and processes that might otherwise remain 
only conjectural | 187| . 



Numerical spacetimes containing a rotating black hole surrounded by an additional sta- 
tionary axisymmetric source were constructed by [|181| 1 (a hole with a thin finite equatorial disc) 



and by [ |237|| (a hole with a thick toroid). To mention just one particular point, [181] ended 
with a prolate horizon (stretching along the axis) in certain cases (when the disc was strongly 
counter-rotating with respect to the hole), which was not observed in |237| 1 . The horizons have 
been known to inflate towards the external sources, so it would be an interesting consequence 
of the interplay between dragging from the hole and from the disc if it were confirmed that they 
can indeed become prolate, even though under extreme circumstances only. 

Note, however, that the main stream of numerical relativity is focused on non-stationary 
problems important in generation of gravitational waves, in particular on black hole inspiral 



The prologue to a major reference |69| (not speaking about its entire content) however convinces us that black 



holes are "the simplest objects" (also "the most perfect macroscopic objects" ) in the universe. Cf. also [ 298 , 115 



3 Within general relativity, there is only one more ingredient that could make the situation even more dramatic: 
non-stationarity. The result would be a gravitational collapse or a collision of already collapsed objects. These 
events are currently devoted an exceptional attention as the strongest sources of gravitational waves. 

4 In the following, we do not review super pos itions, obtained by any of the m etho ds, of black holes with external 



electromagnetic fields. We refer to [p^ p5 |, L9£ , |36| for test EM fields and to [ IOC , |12|, |ll| for exact solutions. A 
more recent account can be found in 
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collision (e.g. [296, p95| , ^] and a thorough review |187f |). 



1.2 Perturbative solution 

A great deal of literature was devoted to perturbations of black-hole spacetimes and several 
formalisms were developed, none of which can be reviewed herej^j In |318|| , a fully explicit 
solution for the perturbation of the Schwarzschild metric by a (rotating) axisymmetric weakly 
gravitating thin equatorial ring was found by solving the perturbed Einstein equations. This 
direct approach is however not practicable for a rotating hole and/or for an extended external 
source with pressure (not mentioning more complicated situations). 

The rotating case (specifically, the algebraically special vacuum case) was seized notably 



by Teukolsky [294] who succeeded in separating the decoupled equations for the first order 
perturbations of a Kerr black hole into the second-order ordinary differential equations for the 
scalar field (for scalar perturbations) or for the Newman-Penrose scalars constructed from the 
electromagnetic field tensor (for electromagnetic perturbations), from the neutrino spinor (for 
neutrino perturbations) or from the Weyl tensor (for gravitational perturbations). By solving 
the "gravitational" Teukolsky equation, the perturbative deformation of the Kerr horizon was 
calculated by [94]. The stationary axisymmetric Green's function of the Teukolsky equation was 
provided by [195]. The most comprehensive expositions of the first-order black hole perturbation 
theory were given by Chandrasekhar |68|, S9| in Newman-Penrose formalism. More recently, [183, 
|214| , [213|| solved the Teukolsky equation in terms of Coulomb wave functions and hyper geometric 
functions. 

In a related method of handling the first-order perturbations of rotating fields, first devised 
for electromagnetic perturbations by Cohen & Kegeles and then generalized to gravitational 
ones [76 1, the perturbation components are expressed in terms of a single (Debye) potential 
that obeys a wavelike equation. Wald showed that this equation is just the adjoint of the 
Teukolsky master equation, while [|160|[ gave a covariant extension of the method of Debye 
superpotentials to neutrino, electromagnetic and gravitational perturbations of all algebraically 
special spacetimes (see [ 30C| , 303] for summary, references and generalization). In the meantime, 
1 77] treated the problem of perturbation of the Kerr horizon after having learnt |76| how to fix the 
perturbed metric from the Newman-Penrose functions. General perturbations of Schwarzschild 



solution were treated in this manner by [301]. Latest advance of the approach consists mainly 
in formulating the conservation laws for perturbations [303, pOfj . 

Let us mention yet another approach [269], also dealing with curvature components but 
not restricted to algebraically special backgrounds. It has been shown, in the static case and to 
the linear order so far, that it can be generalized naturally to self-gravitating matter fields; in 
a spherically symmetric case, the connection with older, metric formalisms of Regge & Wheeler 
and Zerilli was elucidated. 

In the charged case, the perturbation problem introduces coupling between gravitational 
and electromagnetic quantities, implying e.g. conversion between gravitational and electromag- 
netic waves. Interacting perturbations of the Reissner-Nordstrom black hole were studied by 
f32[| (using the earlier formalisms of Regge & Wheeler, Zerilli and Moncrief ) and by |302(| (using 
the method of scalar potentials). A gauge invariant derivation of the basic equations of different 
formalisms was provided in |110| , while [111] related their rotating analog to Teukolsky equation. 



J In monograph [|59[, the author admits that ". . .the account, in large parts, is hardly more than an outline" 
at the beginning of chapter The gravitational perturbations of the Kerr black hole which has 101 pages. 
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The second-order perturbation theory has also been under development since 1970s. For 
Schwarzschild black holes it is surveyed in 



IP- 25(| , while the rotating case is tackled in 57] . 
It should be remarked that solutions describing stationary sources around black holes 
were not the primary aim of the black-hole perturbation strategy. Historically, attention has 
rather been devoted to the stability of black- hole solutions [255|, to the relaxation of black holes 



into a stationary "no hair" state [134], involving the issue of back reaction of a perturbation 

EH) 



IH propagating out of |2(]] or into the hole |24C| , |241| (cf. |l44fl ), to non-stationary processes of 
astrophysical significance such as the behaviour of weakly gravitating particles or waves in black- 
hole backgrounds [253, 148], and, of course, to gravitational waves [231] - presently mainly to 

, |59|, 161]; for a survey, see |296|] ) . 



black-hole collisions as their most prominent source (e.g. 
The latter is exactly an example of a problem where numerical relativity and perturbative 
analysis may 



and should — coexist while mutually checking each other [270, 18, [L9 



1.3 Exact solution 

Perturbative approach is adequate in situations when the external source has only a very small 
effect on the field of the main body. Whenever this is not the case, the superpositions have to 
be described by solutions of the full, non-linear theory. At the end of Jiff Bicak's recent survey 
[33], one is encouraged "not to cease in embarking upon journeys for finding them, and perhaps 
even more importantly, for revealing new roles of solutions already known", because "Is there 
another so explicit way of how to learn more about the rich possibilities embodied in Einstein's 



field equations?" For a more technical demonstration, see the canonical monograph |174| (an 
updated version is awaited); more recent reviews of vacuum fields can be found in [^, |48| . 

The only stationary (in fact static) equilibrium configuration containing more than one 
black hole is the Majumdar-Papapetrou case with extreme centres of the Reissner-Nordstrom 
j, (TSi), where the electrostatic repulsion exactly compensates the gravita- 



type (e.g. g23|, |5C], | 

tional attraction. With rotating sources, gravitomagnetic interaction between spins (repulsive in 
the parallel case — see e.g. [(252fl ) and magnetic interaction between magnetic dipole moments 
(repulsive in the antiparallel case) are also present. Analysis of the momentarily stationary and 
axisymmetric system of two identical sources was carried out by ]96|, [210|| for the Kerr com- 
ponents (with mass M and specific rotational angular momentum a) and by |3j 



211| for the 



Kerr-Newman components (with M, a and charge Q). Equilibrium was found to be only possi- 
ble for the extreme values of the charges {Q = M). In all other cases, supporting singularities 
("struts") occur indicating that a given system of sources cannot remain stationary (or static) 
according to the field equations. Hence, rotating centres can only remain in equilibrium in a 
super-extreme case of two naked sing ularities (that have Q 2 + a 2 > M 2 ) both the magnetic 
dipole-dipole and the gravitational spin-spin interactions being too weak to keep apart black 
holes (having Q 2 + a 2 < M 2 ). (On the other hand, a non-singular equilibrium is possible with 
two spinning Curzon particles |95(| . Most recent discussions of both cases, in particular of the 
character of supporting singularities, can be found in 1 191 , 46, 47|.) 

Accretion discs of astrophysical interest are likely to lie in the equatorial plane of the centre 
and unlikely to bear a considerable charge, so that neither the spin-spin nor the electromagnetic 
repulsion can act within them. One must rather refer to a centrifugal force resulting from the 
orbital motion of the material or to hoop stresses when interpreting the situation. 

It turned out to be difficult to superpose a Kerr centre with an additional axisymmetric 
ring, disc or torus and no explicit exact solutions describing the systems of this kind have been 
found until now. Nevertheless, many stationary axisymmetric solutions of the (electro)vacuum 
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Einstein equations are known that do generalize those containing only isolated black holes. 
These were mostly obtained by indirect methods known as "generating techniques". Two ma- 
jor approaches, developed by the end of the 1970's — the group-theoretic techniques and the 
soliton-theoretic (or inverse-scattering) techniques, work for spacetimes with two commuting 
symmetries. (See [170] for a review and [ |33"| , 84, 



135, 189, 206] for more detailed analysis 
and interrelations between different formulations, e.g. those of Kinnersley & Chitre; Maison; 
Belinskii & Zakharov; Harrison; Hoenselaers, Kinnersley & Xanthopoulos; Hauser & Ernst; 
Neugebauer; Kramer & Neugebauer; or Alekseev.) Other related methods have been proposed 
more recently, e.g. the simplification of the Hauser-Ernst integral equation by Sibgatullin [281], 
"monodromy data transform" by Alekseev |T(| (also |165| ), the static gravitational multipoles 
fl2§| l (cf. 1 2601 ) an d their superposition with stationary fields ] |129| by Gutsunaev & Manko, 
the linear transformation by Quevedo p60f| or the "finite-gap" (algebraic-geometric) solutions 
by Korotkin & Matveev ([172] and references therein). Most of them are briefly compared in 
(appendix 6. of) 111]]. (For other approaches, see e.g. p90|, |82H or the results of Nakamura and 



Kyriakopoulos, referred to and worked out by |293J; cf. also |308[.) 

The crucial point of the soliton generating techniques is a solution of two linear differential 
equations (Lax pair) the integrability conditions of which are exactly the Einstein's equations 
(namely the Ernst equation). The linear problem can be tackled in order to generate new 
solutions from the known ones: given some ("seed") metric with two Killing vectors, it yields 
a different metric of the given type (the procedure is often called Backlund transformation 
in analogy with the technique used for the KdV equation). In such a manner, many known 
spacetimes were reproduced, but also broad families of new solutions were provided characterized 
by arbitrarily large sets of free parameters. However, only a very restricted number of them have 
been given a clear physical interpretation. Though several of these results perhaps represent a 
rotating black hole in an "external" gravitational field (e.g. 261, 208 , 209 , |5l| , fn], [7^]), none of 



the latter has yet been specified to the case generated by a concrete realistic body such as ring, 
disc or torus (cf. the last paragraph of section IV in |190(] ). 

Making the "random" Backlund transformation with some metric, it is, however, difficult 
to require specific physical properties of the spacetime being constructed — one rather "takes 
pot luck" and looks what comes out (see [ |260| | for a hint how to overcome this); section |5.1| 
will illustrate what problems may arise. It would be more appropriate to express the physical 
boundary conditions of the Ernst equation in terms of the quantities which appear in the linear 
problem and then solve the latter. This leads to Riemann-Hilbert problems known from the 
theory of completely integrable differential equations [234, 165, 15]. Tackling the stationary 
axisymmetric boundary value problem with the help of the linear system, the fields of two 
physically relevant types of sources have been discussed: that of black holes [21S , |235|] and that 
of finite thin discs ( |219| , 235, 14, 114] and references therein). There is some hope that the 
above methods could also be used to describe superpositions of both (e.g. |235[| ). 



2 Astrophysical relevance 

The system of a black hole with an accretion disc is very important in astrophysical consider- 
ations. In particular, it plays the key role in models of a whole range of active galactic nuclei 
and of some X-ray binaries. 
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2.1 Galactic nuclei and X-ray binaries 



Extraordinary luminosity generated within a small volume (as manifested by its rapid variabil- 
ity), the presence of very hard (X and gamma) radiation, high rotation speeds and broad velocity 
distribution of the material around indicate that there is an interacting ultracompact object in 
the inner regions of these sources [262, 33, 182, |66|.[] The ambient gas, present in the galactic 
nucleus or overflowing from the double-star primary, typically has enough angular momentum 
not to fall directly (almost radially) onto the object; it rather forms a disc and only gradually 
spirals down. The flow is strongly sheared due to a considerable non-homogeneity of the field, 
which enables viscous torques to heat the gas to high temperatures. (See e.g. [^, |159| , ||, 264 ] 
for various aspects of the modern accretion disc theory. In particular, the most realistic grid 
of non-LTE disc models was constructed in the paper series 

m, B HI EH-) In 

some 

cases, the excess of angular momentum drives outflows through the empty regions along the 
rotation axis of the system. Jets of matter-energy emanate from a number of active nuclei, often 
at relativistic speeds, feeding giant outer lobes distinguished in radio surveys. The formation 
mechanism of the jets is studied by several groups using relativistic-MHD numerical simulations 
(e.g. [16S, ^3|)-0 The high degree of jet collimation evidences an important role of magnetic 
fields [154, |44|, |93|, |218| ] and poses questions concerning the interaction of jets with the ambient 
gas and radiation. (See chapter 9 of | 315 | and [197] for a review.) The entire topic of active 



galactic nuclei is covered in [177] (for earlier references, s ee |317| ) while that of X-ray binaries 
in [192]. Recent status of both areas can be learned from fll52|| . 

In spite of the observational variety of active galaxies (blazars, quasars, Seyfert galaxies, 
radiogalaxies, . . . ), it now seems that we simply look at a similar type of source from different 
directions 



|102| . The discovery of jet structures accompanying several X-ray binaries in our 
Galaxy (the literature speaks of 



''microquasars" [228, |282j]) has strengthened the conviction 
that the stellar-size active sources are also powered by the same accretion mechanism. Further 
evidence is provided by observations of compact X-ray sources in several nearby galaxies |207|. 

The suspicion of the presence of ultracompact objects has also been supported by more 
detailed considerations. Certain spectral features (the iron Ko line in particular) of some active 
galactic nuclei (a celebrated example is the Seyfert 1 galaxy MCG-6-30-15) have been interpreted 
as originating as close to the very centre as a few to a few dozens of Schwarzschild radii,0 so 
they may well represent data from the strongest fields ever met (see [107] for a review, and 
e.g. [266, 267, 155, S8| for models and interpretation issues). The iron lines have already been 
discovered in the X-ray spectrum of several Galactic microquasars ( [|224f| and references therein). 

Also, the frequencies of quasi-periodic oscillations, found at a growing number of stellar- 
size sources, have been identified as the values characteristic for the accretion flow close to a 
black hole or a neutron star [|151| , 230| , 256, 97, 86, 286, 287, 225]. This allows us to consider 



models involving orbital frequencies of "hot spots" in the inner parts of an accretion disc (e.g. 



|153| 1), beat of these frequencies with the centre's spin [ 257[| , blobs on tilted orbits affected by 
Lense-Thirring precession |87|, [l6| , |221|[ , action of a secondary, intersecting the disc repeatedly 
while orbiting the centre [|289([ , oscillations of the (thin) disc itself, perturbed off its circular 



6 In fact, the best evidence for such an object was provided by the nuclei of several "normal" galaxies 171] 
(led by our own [ 121 1 ) which are not so over-exposed with matter and radiation and where it is easier to study 
the motion of individual sources. 

7 For the jet in M87, the first known and the best resolved one, simu lations a re already fitted to observations 
at distances from several dozens of Schwarzschild radii from the centre [150, 109]! 

8 Even a distance of only half the Schwarzschild radius has already been announced (!), giving support for a 
rotating black hole. 
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(equatorial) flow (29|, ^TQ, |l57|] , or waves in tilted discs resulting from the Bardeen-Peterson 



effect [113]. (Clearly, the processes can proceed in a symbiosis, for instance the secondary would 
probably actuate hot spots as well as oscillations in the disc, pulling some of its material out 
into a tilted orbit. 



287, 232 1 for discussions of the 



See [|306(| for a survey and |285| , |86 
individual models in the light of new observations. 

Among the sources displaying the presence of an interacting ultracompact body there is 
a class whose luminosity is still very low and fluctuating within a wide range. This seems to 
point towards the existence of a horizon, because a quasi-stationary inflow onto a neutron star 
should yield a relatively strong and steady output, produced mainly by the gas streams colliding 
with the star's surface. With a horizon, the accretion may (perhaps occasionally) switch over 
to an "advection-dominated" regime when the material plunges in without having released any 
significant fraction of its binding energy 1 233, 119|. See (B| for a general model of advection- 
dominated discs, [141, |142| , 179 , 212 . 204] for recent numerical simulations of their structure and 

for a review. 



spectrum, [217, 25] for a study of jet production in such discs, and |22C 

The main goal of precise measurements of sources suspected of core collapse is a determi- 
nation of three independent parameters — mass of the compact centre M, its specific angular 
momentum a, and radius r where a given process happens. Describing the core field by the Kerr 
metric, the above three parameters appear in the formulae employed in models of the observed 
phenomena. These can be inverted if the necessary independent observables (at least three, in 
general) are really measured; see e.g. [276] for a simple example of such an inference, 120 ] 
for main techniques of fixing the central M in galaxies ( [ |205| proposed a new one), and e.g. 
|l|, H|, @, |20(| for other prospects. 

The dividing line between neutron stars and black holes being given by the Landau- 
Oppenheimer-Volkoff mass limit of about M = (2 4- 2.5)M Q (Mq is the solar mass), attention 
has recently been focused on determination of the rotation rate of the more massive cases. There 
is some evidence that highly active sources (in "hard X-state" , "radio-loud" , with strong and fast 



124]. In particular, the 



jets) host rapidly rotating holes and vice versa [|326| , p7j , |310| , |207j , |217| , | 
popular microquasars GRS 1915+105 and GRO J1655-40 favour rapid central rotation, perhaps 
even close to the extreme value a = M which separates black holes from objects without a 
horizon (naked singularities).^] Such a rapid rotation is expected due to theoretical reasons - 
calculations of the effect of disc accretion on the central body [314, g|. On the other hand, a 

]£] The 



hyper-extreme object would probably be slowed down to a = M by accretion |p6| , 
centre (and the disc) could also lose rotational energy to the benefit of the material outflowing 
from the system. The significance of the Blandford-Znajek mechanism of jet fuelling [ 4l], 263| for 



the active galactic nuclei and X-ray binaries is still under discussion [ 258 , 122 , 198 , 185 1 (while 



it is rather generally accepted for gamma-ray bursts 1 186 1 ) . There also exist other extraction 
possibilities (see [ 309 1 for a review), mainly the Punsly-Coroniti variant of Penrose mechanism 
|25Sj ] (contrary to the Blandford-Znajek one, it does not require a magnetic field threading the 
horizon) which was indeed noticed in numerical simulations [169]. Another efficient scenario has 
been proposed recently [194]. 

9 Note that GRS 1915+105 provided the first di rect evidence of the disc-jet interaction [10C]; GRO J1655-40 is 
likely to have originated in a sup ernova ex plos ion [145]. Their central spins however remain uncertain — cf. the 
disjo int ranges estimated in [124 , [B] and [311]. The spin of the central black hole in our Galaxy is discussed in 
[22G]. 

Conversely, cf. a recently proposed gedanken experiment |9^] which appears to show that it is possible, in 
principle, to turn an extreme (Reissner-Nordstrom) black hole into a (Kerr-Newman) naked singulari ty by a radial 
infall of a (strongly bound, extended and electrically neutral) spinning body. In a related example, |l4C| tried to 
overcharge a near-extremal Reissner-Nordstrom hole. 



7 



Black holes are certainly favoured by the "cosmic censorship" hypothesis [246 1 which 



defends predictability by excluding the naked singularities as the outcome of a realistic gravita- 
tional collapse. However, this "question is still very much open" |248] ] (also ]8l|, |133j[), with some 
amount of evidence as well as counter- indications already given (e.g. [[7(], |l|, |265| , |127| , |149| , 131 [ 
and references therein) whose actual significance is debatable (e.g. [^, |3 1 3| ] ) .f"^ 

The experimental task of deciding between black holes and naked singularities thus brings 
us back to the fundamental problems of the theory. It should be emphasized here that Jiff 
Bicak has always taught us to respect the decisive validity of observations, notwithstanding his 
obvious esteem for the penetrating "pure reasoning" (on the back of an envelope). Jin often 
shows romantic affinity to astronomy and direct contact with nature, though majority of his 
results concern mathematical aspects of the theory. I heard him mentioning the need for a 
balance between the "right and left deviation" in one's work. As a star on the stage of our 
Institute, he also alternates a glowing guitar in "Rock around the clock" and a delicate violin 
in Vivaldi. 

2.2 Gravity of accretion discs around black holes 

It is the subject of a standard university course to discuss the features of the Schwarzschild, 
Reissner-Nordstrom, Kerr and Kerr-Newman metrics (see |33| for a review), exactly describ- 
ing isolated stationary black holes (or naked singularities). Black holes expected in galactic 
nuclei and in X-ray binaries are not isolated, however: matter is concentrated rather than 
absent there. Most of the evidence for ultracompact bodies is in fact based on their interac- 
tion with the surroundings. In theoretical models, the external matter is supposed to form an 
accretion disc around the equatorial plane of the centre (this is justified by the approximate 
reflectional symmetry of the host systems and, at least in the discs' inner parts, by Bardeen- 
Peterson effect |23|; for later results, see |113j |). Accretion is a problem governed by complicated 
magneto-hydrodynamics and radiative transfer which is usually tackled under the assumptions 
of smoothness, stationarity, axisymmetry and simple parametrization of the flow. Gravity of the 
disc itself is neglected, the field is thus fully determined by the centre, either described exactly by 
the Kerr metric, approximated by the Schwarzschild one, simulated by some pseudo-Newtonian 
potential, or quite reduced to a Newtonian (—M/r). 

The main topic of the present paper is an inclusion of the disc gravitational field in a fully 
general relativistic description. Newtonian treatments of the problem began with Ostriker's equi- 
libria of uniformly rotating, polytropic, self-gravitating slender rings [242]. The paper opened 



with references to previous centuries and concluded with remarks on extragalactic radio sources 
and stability of rings, both of which were to boom soon. Later, polytropic self-gravitating con- 
figurations with a realistic equation of state and opacity were calculated by |173|1 (cold case) and 
|11£, 42, |132 | (hot case), already suggesting the possibility of a black hole present in galactic 
centres. |227 ] incorporated viscosity (in a self-similar case).[^] In the meantime, basic formal- 
ism for studying self-gravitating discs around black holes in a stationary axisymmetric general 

11 Black holes are sometimes said to reside where God divided by zero. Regardless of whether he or she used 
Good coordinates or not, this is consistent with his/her b eing Almighty. Certain current readings of the cosmic 



censorship conjecture ("God abhors a naked singularity" [135]) sound somewhat disturbing, on the contrary. 

12 The papers must be studied in detail in order to understand precisely what was done. The topic of accretion 
discs is very wide and conclusions found in the literature are often hardly comparable due to different kinds 
of simplifications made, relating not only to the gravitational field and symmetries, but mainly to the rotation 
law, accretion rate, equation of state, viscosity, heat production, radiation pressure, opacity, boundary conditions 
(material supply), etc. 
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relativistic case was given by p2], 112 



Though the mass of real inner accretion discs is usually claimed not to exceed few percent 
of the central mass,P| the role of their self-gravity is not fully understood yet. Are any properties 
of accretion discs so sensitive to the details of the field that even a small effect of the disc itself 
can alter the flow significantly? It was shown by [ 28C ] that the very global structure of the disc is 
an example of such properties. This paper presented an extension of the disc "a-models" (both 
in their Newtonian version of Shakura & Sunyaev and in the relativistic version of Novikov & 
Thorne) which includes the effects of self-gravitation, radiation pressure and variable opacity 
source. The self-consistent configurations thus obtained were considerably thinner than their 
test counterparts, showing sharp radial decrease of the vertical dimension in the middle region. 
The implications of self-gravity concerning global properties of the disc (e.g. distribution of 
angular momentum and outer radius) were discussed and interpreted in a Newtonian analysis 
Q, the roles of the disc density and total mass were distinguished and specified in particular. 
With increasing gravity of the disc, regions occurred where free circular motion became unstable 
or even impossible, thus it was stated that "astrophysically relevant models of thick accretion 
discs must incorporate self-gravity effects in order to be self-consistent" (the same had already 
been indicated in |280| , |316| ] before). 

Changes in the properties of circular equatorial motion were later observed in the pseudo- 
Newtonian study ||162|] of a rotating black hole surrounded by a massive thin equatorial ring, in 
a relativistic analysis of the corresponding static case |34|] , and also in the perturbative solution 
1 315] for a slowly rotating, weakly gravitating thin ring around a Schwarzschild black hole 
1 31£ 1 . These changes were claimed to result in modifications of the observable characteristics, 
e.g. in higher temperature of the disc. The latter was confirmed by numerically constructing 
equilibrium structures of self-gravitating polytropic thick discs in a pseudo-Newtonian potential 
|42| . The accretion-disc "self-field" was also discussed in connection with the disc's vertical 
structure and stability. It turns out that certain modes of instability can be damped whereas 
other are amplified by the disc's gravity (see [244, |l3|). The latest fully relativistic results - 
employing numerical spacetimes obtained in 1 181 , 237| — include the computation of emission- 
line profiles from self-gravitating discs around rotating black holes; the case of light, thin finite 
discs was treated by [156] while that of heavy toroids by [305]. 

Self-gravitating accretion discs also appeared in other contexts, some of which are inter- 
esting from the point of view of general relativity (e.g. heavy discs around white dwarfs as a 
ponderable stage of a double-white-dwarf merger [|130| ), whereas other are not relevant for the 
present review (young stars [32C, |321fl , outer parts of nuclear discs, galactic discs [27£], effect on 
bipolar jets [ 117] ) . 

Recently we have studied exact solutions for a static (Bach-Weyl) ring |64j and static thin 
annular disc |188| around a non-rotating black hole; some of the results are presented in section 



We also checked the possibility to generalize them to a stationary situation | 325 | (section ||[) 
when the hole and the disc were allowed to rotate and the self-gravity effects could be even more 
pronounced. Actually, in relativity, the kinetic (in fact any) energy also generates the field and 
mass-energy currents give rise to frame dragging. Whereas mass of the hole should be dominant 
in a real accretion system, the disc can bear much (even most) of the angular momentum, thus 
modifying the gravitomagnetic field of the centre significantly. This could be important for the 
mechanisms of rotational energy extraction. 



This only applies to the ultracompact centre (cf. e.g. [316, 231], however), not to protostellar discs 1 321 ] 
where even the opposite may well happen [320]. 
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2.3 Gravitational collapse and gamma-ray bursts 



We should not forget the accretion system where the mass of the disc or torus can actually be 
comparable to that of the central hole. It is considered in the literature as a transient state of 
gravitational collapse that could power the gamma-ray bursts. 

In a present-day relativistic astrophysics, perhaps the most popular system is a close 
binary of ultracompact objects. The famous Hulse- Taylor binary pulsar (PSR 1913+16) turned 
out to spiral in at the rate predicted by general relativistic loss of angular momentum via 
gravitational radiation [291 1, thus providing the first, indirect confirmation of the existence of 
gravitational waves. Nowadays late stages of evolution of such systems (inspiral, merger and 
"ringdown") are studied numerically or perturbatively by a number of research groups, mainly 
as the most promising source of gravitational waves [^96| , and also as an engine suggested 
for the gamma-ray bursts [223]. In particular, at late stages of the neutron star-black hole 
binary inspiral, the neutron star tidally breaks up. Loosely speaking, if this happens outside the 
innermost stable orbit around the hole, the neutron debris spreads into a strongly magnetized 
torus (otherwise it plunges into the hole). Similar outcome is expected in case of a binary with 
two neutron stars (which is expected to be less frequent than the black hole-neutron star case 
|29|] ). It was pointed out in numerical analyses that the black hole-neutron torus system was 
liable to runaway instability (see f216j l and references therein). Depending on the exact shape 
of the created magnetosphere, the energy release is mainly ensured by leptonic winds from 
pair-creation along the axis of rotation or /and axially collimated Poynting outflow resulting 
from the Blandford-Znajek mechanism |18G| , |307[ . In |[307|] , "Long/short gamma-ray bursts are 
identified with suspended/hyper- accretion onto rapidly /slowly rotating black holes. . . .In long 
bursts, the torus is expected to radiate most of the black hole luminosity in gravitational waves. 
This predicts that long gamma-ray bursts are potentially the most powerful burst-sources of 
gravitational waves in the Universe." (This thorough review also conjectures about connections 
between gamma-ray bursts and active galactic nuclei and micro-quasars.) 

The black hole with a disc/torus can also arise as the outcome of a gravitational collapse 
of a massive rotating star. This was already mentioned in [231], while more recently in scenarios 
of "failed supernovae" from the cores of young stars, either isolated |205H or in a binary (this 
case is referred to as "hypernovae" ) p43|, 54]; for a review, see discussions in [307, 



3 Stationary axisymmetric spacetimes 



It is natural first to try to include the gravitational effect of the external source in the simplest 
case when the system is stationary and axisymmetric. Stationary axisymmetric spacetimes are 
of obvious astrophysical importance: they describe the exterior of bodies like stars, galaxies or 
accretion discs in equilibrium. The issue has been exposed in many places (e.g 
74] or chapters 17-19 of 



192, 31 



318 



55 



103, 14 



particular solutions, see |45 



In the Weyl-Lewis-Papapetrou coordinates (t, p 
can be written as0 



74]]); for more information about the most important 
z) of the cylindrical type, the metric 



ds 2 



-e 2u dt 2 + p 2 B 2 e 



1v 



(d0 - uj&t) 2 + e 2X - 



■ 2v {dp 2 + dz 2 ) 



(1) 



14 The quantities will be given in geometrized units in which c = G = 1; the signature of the metric tensor g^ v is 
( — I — | — h). Greek indices run from to 3 and Latin indices (i, j, . . . ) run from 1 to 3; indices from the beginning 
of the Latin alphabet (a, b, ... ) represent cyclic coordinates i and (f). Partial differentiation is denoted by d or 
by a subscript comma, covariant derivative is denoted by V or by a subscript semicolon. 
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where the unknown functions v, B, uj and A only depend on p and z; uj is interpreted as the 
angular velocity of inertial space rotation with respect to observers at rest in spatial infinity. 
The metric coefficients g a b have invariant meaning, they can be expressed in terms of Killing 
fields = dx^ jdi and = dx^ /d(f>: 

gtt(=-e 2u + uj 2 g (f)<f) ) = k i k\ (2) 

g t< t> (= -ujg^) = hm\ (3) 

g H (= p 2 B 2 e~ 2v ) = m.m'. (4) 
The simple result for the corresponding subdeterminant, det(g a b) = —p 2 B 2 , implies 

g tt = -e- 2v , g t + = -e- 2v u>, g** = -e~ 2v u 2 + p- 2 B~ 2 e 2v . (5) 

3.1 Field equations 

The Einstein equations read 

V ■(pVB) = 87r P B(T pp +T zz ), (6) 

V • {p 2 B z e~ iv Vu) = -\^Be 2X - 2v T\, (7) 

V • (BVis) = ^p 2 B 3 e- 4u (Vu) 2 + 47rBe 2A (2T tt + e~ 2v T), (8) 

where T = T° and V and V- stand for the gradient and divergence in a flat three-dimensional 
space with coordinates (p,(f>,z); thus VX = (X >P ,0,X ;Z ) and V • X = p~ 1 [(pX p ) >p + (pX z ) tZ ] in 
the axially symmetric case. Once B, uj and v are known, A can be integrated from equations 

B\ z -B )Z + p (B tP \ z + B jZ \ p - B )PZ - 2Bv p v z ) 

+ ^p 3 B 3 e- 4u u p u z = 8tt P BT pz , (9) 



2BX tP - 2B :P + p 



1 



2-B p A )P - 25 2 A ;2 - B iPp + B t zz - 2B(v 2 p - v 2 z 



+ -p^e- 4 >f p - u 2 z ) = 8tt P B(T pp - T zz ). (10) 

Although a suitable combination of solenoidal motions can also satisfy the assumption of 
axial symmetry, it is natural to expect that the axially symmetric source just rotates in the 
azimuthal direction, with an angular velocity O = d(p/dt. The corresponding four- velocity is 

= u t (k' x + nm t ") =u\l, 0,0,0), (11) 

where (u*) -2 = — gu — 2VLgt<f> — ^g^ = e 2u (l — v 2 ) and v = pBe~ 2v (Sl — uj) is the linear speed 
with respect to the local zero-angular-momentum observer. 

For a source made of ideal fluid, having a total mass-energy density e and pressure P as 
measured in a co-moving system, the energy-momentum tensor is T pv = (e + P)u p u u + Pg pv 
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and the equations of motion read (e + P)a fl = — P M , where the four-acceleration a p = u p - u u u 
can, for instance, be written as 



a. 



2 —^~ - 2e 2v (\-v 2 ) 

_ 2e 2i V M + 2(0 - ujg^u^ - (0 - lo) 2 'g^ lfl 
2e 2v {l-v 2 ) 

From equations @-(|8]) one then obtains 

V • (pVB) = 16tt pBe 2X ~ 2v P, 



V • (p'B'e-^Vu) = - 16^05 V A ~ 4iy (e + P)- 

V • {BVu) = ^p 2 B 3 e- 4u (Vto) 2 + ATiBe 2 ^ 2 " 



( e + P)i±| + 2P 



and A is given by 



B 



(pB)l + p 2 B 2 



\b, p + p b-\b 2 p + b 2 



+P 



p z B 2 e~^ \{pB), p (uj 2 p - lo 2 z ) + 2pB jZ u z u, 



B 



( P B) 2 +p 2 B 2 \ B ' Z + 9 [ pB ' z ( U >* ~ V J + 2 (P B ),P v ,P u r 



- \p i B 2 e~ iu [pB tZ (u? z - u 2 p ) + 2{pB) tf> u, p u z 



(12) 

(13) 
(14) 

(15) 



(16) 



(17) 



The unknown metric functions are subject to boundary conditions on the horizon (if there 
is one), on the symmetry axis and at spatial infinity. The metric must be regular on the horizon 
and on the axis; in the case of an isolated source one requires asymptotic flatness. The conditions 



were discussed by EL 22 j, for example. 



3.2 Vacuum case and the Ernst formulation 

In the following sections, we will be interested in vacuum solutions or in those with infinitely 
thin sources. In both cases P = and the equation ( |l~3| ) has only one solution with a satisfactory 
asymptotic behaviour [namely (pq)], B = 1. Only two of the field equations then remain, ([f],||), 
in the form 



V • (p 2 e- iv Vu>) = 0, 

V72 1 2 -Av/^t \2 

V v = -p e (Vw) , 



(18) 
(19) 
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and the relations (|ll|[L7]) reduce to 

1 

,p ",zj V~,p 



\ P = m - v\) - -y^^i - 4), (20) 



1 

\ z = 1pv^v^ z - -p 3 e~ u u >p oj >z . (21) 

Many other formulations of the stationary axisymmetric problem are possible, starting 
from a different definition of the four (or three) unknown metric functions. In the vacuum case 



the formulation by Ernst [104, 105] is often considered, where the field equations are translated 
into the (Ernst) equation 

- g tt V 2 £ = (V£) 2 (22) 

for a complex (Ernst) potential £ = —gtt + whose imaginary part ip is given by 

pip,z = gtt9t<t>,p ~ 9tt,p9t<t>, pip, P = 9u,z9t<t> ~ 9tt9t<t>,z\ (23) 
A is found by a line integral as above. 

3.3 Horizon 



(See e.g. [31, 6^, 115] for thorough accounts.) A black-hole horizon is a null hypersurface below 



which the spacetime is dynamical. In Weyl-Lewis-Papapetrou coordinates, it is located where 
det(g a b) = —p 2 B 2 is zero. Also vanishing there is the lapse (or redshift factor) a = e v or the 
magnitude of the Killing bivector k\ K mx\; this is clear from relations 

- 2k [K m x] k^m x] = - det(g ab ) = a 2 g Hl (24) 

because g^ (as well as g pp ) must be regular (and positive) on the horizon. In the vacuum case 
B = 1, so the horizon lies on the axis (p = 0). The interior of the black hole thus has to 
be studied in different coordinates, e.g. in spheroidal coordinates of the Boyer-Lindquist type 
(t,r,6,4>), introduced by the transformation 



p=VA sin 9, z = (r- M)y, (25) 

where A = (r — M) 2 — k 2 and y = cos#, M being a scale parameter (it usually represents the 
total mass) and (±)fe determining where the horizon reaches up along the z axis. Sometimes an 
isotropic radial coordinate R is employed, given by 

A = fl2 ( 1 -^) 2 - < 26 > 

On the horizon, A = and both radial coordinates are constant, r = M + k = rn (hence, 
z = ky), R = k/2 = Rr (there are more black-hole horizons in general, we mean the outermost 
one here). 

There are three important parameters of the horizon: its surface area (proportional to the 
entropy of the black hole), the surface gravity (proportional to the horizon temperature) and 
the angular velocity relative to infinity, 

A = 2ir J \J (geeg^u d6>, k = [Ve v | H , w H = oj(r = r H ); (27) 
o 
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k and wh are constant all over the horizon. 

For spacetimes with a symmetry lower than spherical, the question of the "true shape" 
of the black hole, not distorted by coordinates, is non-trivial. How is the hole influenced by 
rotation, electromagnetic field or an external source? The answer is obtained by representing 
the horizon as a two-dimensional surface in a three-dimensional Euclidean space. Following 



1 282], the two-dimensional metric is first rewritten as 



&s 2 = ^[h- l {y)dy 2 + h{y)dct>\ (28) 

where h(y) = (4ir/A)(g ( j )( f ) )i{. An isometric embedding of the two-surface (y,(fr) in IE 3 with 
coordinates (X, Y, Z) is then given by 

y I 

X = ^Vhco S cj>, y = A^sin0, Z = ±jJ\(l-\h£)ty. (29) 

o v 

If its Gaussian curvature 

Cn = --^-hyy (30) 

is negative anywhere, the horizon cannot be (globally) embedded in IE 3 . It was shown in |283| 
that rotation can flatten the horizon considerably, the Gaussian curvature of a Kerr black hole 
finally turns negative at the axis for a/M > \f?>j1 (a denotes the angular momentum per unit 
M). Some of the black holes swirling in galactic nuclei may even be impossible to imagine! 

3.4 Static case 

In the static case lo = and (|l]) acquires the Weyl canonical form 

ds 2 = _ e 2, dt 2 + p 2 e -2u d(/) 2 + e 2A-2, (dp 2 + ^2) (31) 

More precisely, the metric can be written in this way provided that + T* = (which is 
fulfilled by incoherent dust, certain electromagnetic fields or certain infinitely thin sources, but 
not by a fluid with non-zero pressure), otherwise there remains B 2 in g^. 
Equations (|6|)-([l0|) reduce to the Poisson equation 

V 2 v = p~\ p + u pp + u zz = 4vre 2A -^(T* - T*) (32) 
for v and to relations 

\z ~ 2pv, P v,z = 8irpT pz , A iP - p(v 2 p - v 2 z ) = Airp(T pp - T zz ) (33) 

for A. 

3.5 The axis and the equatorial plane 

On the symmetry axis (p = 0), the regularity of metric requires e A = B, so 

&s 2 = -e 2v &t 2 + B 2 e- 2 »dz 2 . (34) 

The formulas can also be expected to simplify in the equatorial plane {z = 0) if the spacetime is 
reflectionally symmetric with respect to it. This is usually the case in astrophysically motivated 
considerations. As a special case, the solutions with stronger, cylindrical symmetry (where 
another Killing field, dx^/dz, exists) can be mentioned, where any of the planes z = const is 
equatorial in this sense; however, there is no place for black holes in these spacetimes. 
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3.6 Singularities 



From astrophysical point of view, a given solution is problematic mainly if it contains physical 
singularities on or above the horizon, if it shows bad asymptotic properties or if it does not 
correspond to any realistic source. The existence of a singularity follows — according to its type 



(see [101, |80|) — e.g. from the divergence of scalars constructed from the metric tensor and 
its derivatives or from the divergence of physical (tetrad) components of the Riemann tensor. 
(Usually, the Kretschmann invariant R^ up(T R^ vpa is checked first.) However, the spacetimes are 
known which do contain singularities, although all their curvature invariants vanish. It is thus 
difficult to verify the non-existence of singularities; the only proof of regularity at a given point 
is to find a (local) coordinate system in which the metric is smooth enough there. 



3.7 Analytic extension and global structure 

Astrophysically relevant is that part of spacetime which can communicate with spatially remote 
regions, i.e. the "region of outer communications" outside the horizon. General relativity is 
interested in the black hole interior as well — in fact even there an observer can perform physical 
observations (and perhaps even survive). Unfortunately, in Weyl-Lewis-Papapetrou coordinates 
the metric is not defined below the horizon (this would correspond to imaginary radius p) and 
one has to extend it there. Such a task may be solvable but numerically and hardly expectable 
regions may open. What is the type and shape of the singularity under the horizon? How does it 
change when an additional source is "switched on" above the horizon? The answer is unknown 
even for quite simple superpositions (on the other hand, in a highly symmetric case it may well 
be possible to solve the dynamical problem |178| j). 

Far away from the horizon, one inquires for the asymptotic behaviour of the field. That 
of an isolated stationary source falls to zero there in a specific manner — the spacetime is said 
to be asymptotically flat (see e.g. 1 312 1 , chapter 11, or |26|], section 3). In case of the metric ((!]), 
it must hold, for r — > oo, 

total mass _ . _ 2 . .„_. 

v = + 0(r 2 ), (35) 

r 

B = l + 0(r- 2 ), (36) 

2 • (total angular momentum) _ 4 
u= -3 + 0(r ), (37) 

A = 0(r" 2 ). (38) 

The courses can be expressed, in the same manner, in terms of the isotropic radial coordinate R 
or in terms of the radius \J p 2 + z 2 , since R = r — M + 0(r -1 )^] and \J p 2 + z 2 = r — M + 0(r~ 1 ) 
for r — > oo. 



3.8 Circular orbits and ergosphere 

The simplest type of motion in stationary axisymmetric spacetimes is that along spatially circular 
orbits {p = const, z = const) with steady angular velocity £1: such a motion takes place 
along the symmetries, namely tangent to each individual circular orbit is a Killing field (one 

15 This is clear from an explicit form of (||), R = |(r - M + \/A), or r = R + M + 
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speaks of "quasi-Killing" trajectories); physically speaking, an observer on a circular orbit sees 
time-independent field around and is thus called the stationary observer. In order that the 
corresponding four-velocity (O) lie inside the light cone, fl must fall within the values 



OL 

Omax = u) db — ■ ; (39) 



mm 



on the horizon (where a = 0) the permitted range narrows down to a single value wh- 

Calculations involving acceleration are more difficult as they also contain derivatives of 
the metric. The circular-orbit four-acceleration (|l2| ) has at most two non-zero components (a p 
and a z ); /x-component is zero if 



si = ll± = , (40 J 



where the upper/lower sign corresponds to a "co-rotating" /"counter-rotating"^] orbit. 

In the equatorial plane, a z = due to the symmetry. If = £l± (here with p = p), a p = 
as well and the worldline is a geodesic. Among equatorial geodesies, three particular cases are 
most important: the photon, the marginally bound and the marginally stable orbits. The photon 
orbits bound the regions where circular motion is time-like; they are given by equalities 

fi-fc = ilrnax (41) 
min 

[p = p being the only non-trivial index value in ( pjO| ) now] . The marginally bound orbits limit the 
regions where particles on circular orbits have lower energy than is necessary for the existence 
at spatial infinity; the limiting case is given by 7(fi = Q±) = 1, where 

7 = -u t = -vtfett + Slguj)) =u t e 2u + ul and I = = u*s^(il - u) (42) 

are specific energy and specific azimuthal angular momentum at infinity. The marginally stable 
orbits bound the regions where circular motion is stable with respect to perturbations within 
the orbital plane; they satisfy [£(p, = Q±)] tP = 0. They play a key role in the theory of accretion 
discs: the matter is "swept away" from unstable sectors; in particular, the innermost marginally 
stable orbit should represent the inner disc rim. (We will see in section || that perturbations 
perpendicular to the orbital plane can also be important.) 

A peculiar feature of rotating fields are the dragging effects. Their history starts from 
the famous Newton's bucket experiment (if not earlier, see Dolezel in this volume) and contin- 
ues notably by Mach's ideas of relativity of motion and inertia towards Einstein and Lense & 
Thirring who showed, in the early times of general relativity, that the effect is present in this 
theory (see the 1913's letter of Einstein to Mach ]166| ).P1 The effect has an obvious analogy in 
electromagnetism, where (electric) currents generate the magnetic component of the field. This 
relationship was already noticed by Einstein (see his Prague paper of 1912 fl67f| ) and today the 
gravito-electro-magnetic approach to rotating fields prevails. 

An extreme implication of dragging is the occurrence of the ergosphere in the vicinity of 
ultracompact rotating objects. In this region r2 m ; n > 0, thus the stationary observer cannot 
be static (at rest relative to infinity), albeit he still withstands the radial attraction. The 
static-limit surface, given by O m i n = and hence by gu = 0, also represents the set of points 

16 These terms may be misleading in spacetimes with large angular momentum. 

17 Mac hian i nspiratio n keeps vivid within contemporary relativity — see, for example, the last-decade references 



[ 14^ , |7|, [Tel or the whole conference Q. 
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from where signals emitted by static observers (0 = 0) reach infinity with an infinite redshift. 
Ergospheres received considerable attention (and the name from epjois) after processes had been 
suggested (working only there) [ |246| ] by means of which the rotational energy of the centre could 
be extracted (without violating the second law of black-hole dynamics that forbids the horizon 
area to decrease). 

With the metric ([[]), gu = corresponds to pBuj = a 2 which can only be satisfied by 
p > 0, so the static limit really lies outside the horizon;^ in particular, p = implies a 2 = 0, so 
the static limit touches the horizon at the axis. 



3.9 Sources 

In order to complete the relativistic solution, one must also describe their interior (where T^ u ^ 
0). Finding a realistic interior solution is a remarkable achievement in general, mainly if it 
should match smoothly the vacuum exterior (e.g. [18C, 24£, 73] and references therein); this 



is also confirmed by the history of searching for a source of the Kerr field (see [{38| , |236| and 
references therein). 

The problem is much simpler if the source is infinitely thin: the solution is then vacuum- 

— 1/2 

type everywhere, with the energy-momentum tensor T^ v = g Z z S^ v 5{z) found from the jump 
of the normal field across the source like in electrodynamics; the appropriate covariant method 
is known as the Israel's formalism |24|] . It starts from a projection of the exterior (here vacuum) 
metric g^ v onto a three-dimensional hypersurface S, representing the history of the source: the 
projection from one (+) and the other (— ) side of S yields the induced three-metrics 

± 9 A B = 9„ ± e^e^ (43) 

where < e*t \ are bases tangent to S. The hyperplanes S are however the back and 

I AJa=0,1,2 to Jr r 

the face of the same surface, so we can choose = ~e^~ and only the normals and ~r*/ J 

' A A 

differ. The surface density of the energy-momentum tensor Si a reads 



SAB = ^r\\ K AB\-9AB\ K( 6\)' (44) 

where [-K^] denotes the jump of the exterior curvature across S, i.e. [-f^g] = + ^ab ~ ^-Ab^ 
with 

±K AB = CA eU B V » ±n » = -^^V^. (45) 
The case of a thin source (a disc) in a stationary axisymmetric spacetime has notably been 



studied by [184, 126]- In a disc lying at z = const the radial pressure is zero, S pp = 0, and for 
a disc in the equatorial plane [z = 0) of the renectionally symmetric spacetime S zz = 0, too. 
Only Sfrh is non-zero off the diagonal. In the Weyl-Lewis-Papapetrou coordinates the formula 
dH) can be written as ( fll&l , chapters 3 and 4ff] 

S ab = -4^(—) ■ (46) 




18 In coordinates which also describe the black hole interior, two horizons and two static limits are found in 
general. We mean the outer horizon and the outer static limit everywhere. 



19 The indices a, b, ... only run through t and <fi as anywhere. 
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Using (|5|) and (21), this yields 



si 



57T 



4^(1 - pv tP ) - p e v u )Z {uj - poJ iP ) 



8tt 



-p 2 e^u 



-x 



8tt 



Aujv 



i , 2 -4z/ 2n 

1 + p e u )u) j 



8tt 



4z/ p ^, z - pe "o^fo; + pto 



(47) 
(48) 
(49) 
(50) 



The total mass and angular momentum of the disc are fixed by the Komar integrals that 
lead to 



M = l[ 9 * 9tc , z pdp = -\[ 9 ^(^) dp, 
2 J 4 7 p \gwj 



b 

oo 2 



J = \ J 9 tC 9c<t>, z pdp =\J —u, z dp 



(51) 



(52) 



where p = b > is the position of the disc inner rim. All the z-derivatives are understood to be 
calculated in the limit z — > + . 

In the event that there is a black hole present in the centre of the disc, its contributions 
Mn and Jh need also to be included; for a given type of solution, the aggregate parameters of 
the spacetime can be written as the sums Mh + M. and Jh + J (see e.g. |62|]). The black- hole 
mass is related through the Smarr formula (162], equation (6.297)) 



Mh = 2wh Jh + -r- = 2^h Jh + k 
with the angular momentum 

7T 

J H = -i |(A 2 e - 4 ^, r ) H sin 3 OdO. 
o 

The last step is a physical interpretation of the source. First, if 

2v-2X 



16tH 



is not negative, then there exists a stationary observer with the angular velocity 

2e iu v : 



n lso = (2^)~ x (s* - St - Vd 



oj 



p 2 u t 



- + 



\\7 



■oj 



(53) 



(54) 



(55) 



(56) 



with respect to whom the energy- momentum tensor assumes a diagonal ("isotropic") form 
S» v = umf so < so + P^ < , where uf so = nf so (l, 0, ft iso , 0) and v? BO = -p _1 (A so , 0, 7 iso , 0) are 
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the observer's four-velocity and unit base vector in the ^-direction, £j so and 7i so being the cor- 
responding specific angular momentum and energy at infinity (see section |3.§D-P| The observer 
measures the surface density 

2 att _ p2 M 

ut h- +n- I ) [ ' 

""isov fiso \ ^'isomso^ 

and tangential pressure 



/ISO \ J 'lso-wso 

If w > P > 0, the energy-momentum tensor can represent two equal streams of particles moving 



on (accelerated) circular orbits in opposite directions at the same speed yP/w. This "speed 
of sound" was shown in |126| to be equal to the geometric mean of the local prograde and 
retrograde circular geodesic speeds v± as measured by the observer u^ so : 

\v+V-\=P/w, (59) 

where 

(fiso)t + {Viso)^± 

v ± = 7 s 7 7n (6°) 

(Miso)i + («iso)</.i2± 



and Q± are the angular velocities of the prograde and retrograde circular geodesies (|40| ) . 



It was shown in [184] that in the absence of radial pressure the source disc can be com- 
posed of two counter-rotating streams of non-interacting particles on time-like circular geodesies, 
provided that the following conditions are satisfied: both geodesic frequencies (|40|) lie within the 
sub-luminal interval (r2 m m,Q max ), the tensor fl46]) yields S a bX a X b > for at least one two- vector 
X a , and det(S a b) > 0. The last two assumptions ensure that the superposition 

S ab = w+u a + u b + + w^u a _u b _ (61) 

of the two dust components on circular geodesies u± = it±(l, fi±) comes out with positive proper 
surface densities of the streams 

** = ± MV + -«-V (62) 

In a static case, the above formulae simplify to f2i so = 0, 



„u— A ,,v— A 



W 



s t = (1 - PKp) v ,*i P = S t = ^—P",pV,z, (63) 



~2 2vr v ' ' F/ f 2vr 



DO 



M = / i/ *pdp, .7 = 0, Mh = k = M, J h = 0, (64) 



6 



«± = pe" 2 ^ ± = ±J-^~ = ±y^P/w . (65) 



2(1 



The case of D < (not involved here) represents discs with non-zero heat flow. Then the energy-momentum 
tensor must be written in a more general form S 1 *" = wu^u^+Pv^v^+K^^v^+v^u^) with K = yJ—D/2, 
and n iso = {2S t 4> )- 1 {S^ - Sf) = u - '^f- See p 
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Hence, the disc particles circumscribe circular geodesies in the spacetime which they generate 
themselves. 

Now further physical requirements can be raised, namely the weak energy condition (w > 
0, P > —w), the dominant energy condition (w > 0, \P\ < w), and the non- negativity of pressure 
(P > 0). Their combination yields w > P > 0. The requirements are very simple in the static 
case: if v z {z = + ) > (which is fulfilled for a realistic matter), then P > is equivalent to 
v ,p > which is satisfied below the Lagrangian point of zero field (where u„ = and v± = 0). 
In regions with tension {P < 0) , the particles are more attracted by the outer parts of the disc 
than by the central black hole, thus no circular geodesies exist (then hoop stresses have to be 
employed to interpret the disc matter); this typically happens in discs with too much matter on 
larger radii. The other condition w > P is equivalent to v± < 1. Hence, the above constraints 
are in fact ensured by the obvious condition < v± < 1 which can be written explicitly as 
1 > > 1/2. 



4 Static thin discs around non-rotating black holes 

In a static spacetime rotation is excluded or it has to be compensated exactly as in the case 
of counter-rotating streams of matter. Within the vacuum outside of the sources, the static 
axisymmetric Einstein equations (^) and (|33| ) yield the Laplace equation for v and a simple 
quadrature for A (calculated along a vacuum path going from the axis to the given point), 

p,z 

V 2 i/ = 0, A= J p[{y 2 p - vl)dp + 2v p v z dz}. (66) 

axis 

Thanks to the linearity of the Laplace equation, the solutions can simply be added to obtain 
fields of multiple sources. However, such a superposition may not be physically acceptable 
- supporting singularities ("struts") may occur in calculating A or it may not be possible to 
interpret the resulting system of sources in a reasonable way (for instance, the matter works 
out with negative density or pressure or moving at a superluminal speed). Consequently, even 
the static axisymmetric case has only afforded a few realistic superpositions. Let us refer to two 
of them involving a Schwarzschild black hole: the superposition with an infinite annular thin 



disc, obtained by Lemos & Letelier [185] by inversion of the first counter-rotating finite disc of 



Morgan & Morgan [229], and the one with an inverted isochrone thin disc, obtained by Klein 



1 164]. The Lemos-Letelier solution has turned out to allow for physically satisfactory situations. 



We briefly summarize its properties below. 



4.1 Superposition with the inverted first Morgan-Morgan disc 

The inverted first Morgan-Morgan counter-rotating disc has a surface density 

<p) = 2 ^4^^ ( 67 ) 

Ai and b denoting the mass and the Weyl inner radius of the disc. The matter is thus concen- 
trated near the rim (w is maximal at p = 2b/y/3 = 1.1156) which is supposed to be the case 
in real accretion discs. The resulting spacetime is given by a sum of potentials which in Weyl 
coordinates read 

1 d!+d 2 - 2M 

^ = 2 ln d 1+ d 2 + 2M (68) 
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and 



M 



vr ( p 2 + z 2y/2 



2p l + 2z 2 



> p 2 -2z 2 
p 2 + z 2 



arccot \ 



lo-{ P 2 -b 2 + z*) 
2(p 2 + z 2 ) 



to 0,2^ 2^ 2^ a-(p 2 -b 2 + z 2 ) 
- {3a - 3b +P +Z) V 8(p 2 + , 2 ) 



(69) 



where d\ 2 = \/p 2 + (2 =F M) 2 and <r = \/ (p 2 — 6 2 + z 2 ) 2 + 4b 2 z 2 ; M is the Schwarzschild mass. 
On the axis (p = 0), fdi sc reduces to 



2.M 



7T Z ■ 



[z 2 + 6 2 ) arctan — — b\z\ 



while in the equatorial plane (z = 0) 



^disc(P > b) 



M 



2p 2 



(70) 



(71) 



^disclP < 0) = 

7rp 



' 6 2 \ p 6 2 

2-^larcsm- + W^-l 



(72) 



In [277] we illustrated the shape of the superposed field by drawing its field-lines, namely 
the integral curves of the four-acceleration = V M z/ of a static congruence. Then we plotted 
the flattening of the black hole with increasing relative mass of the disc or with decreasing inner 
disc radius. The explicit (static axisymmetric) form of the horizon Gaussian curvature (|30|), 



Ch 



1 - 



ext,j/ 



:i 



y 2 )(2 K 2 



cxt,y 



u ext,yy) 



4_/\^2 e 2!/ cxt (j/)-4^ cx t (1) 



(73) 



(the exterior potential u eX ^ is represented by fdi SC in our case), implies that starting from certain 
(though unrealistic) values, Ch becomes negative at the axis. 



The next paper [275] dealt with the motion of free test particles in the superposed fields, 
especially with the influence of the disc parameters on the positions of important equatorial 
circular geodesies (the photon, the marginally bound and the marginally stable one) which are 
critical for the disc accretion flows. The following questions arose: with increasing disc mass, 
how does the inner rim have to shift if the whole disc is to be stable permanently? Does it shift 
towards or away from the horizon starting from the Schwarzschild radius r = 6M where the 
marginally stable orbit lies in the case of a test disc (A4 = 0)? How do other properties of the 
resulting spacetime evolve? The answers were given in papers |274| , |275 |, |32z 



A sequence of superposed solutions was generated (figure 4 in [|274j ] ) by gradually increasing 
the relative disc mass A4 /M, while demanding that (i) all the disc matter can be interpreted as 
two equal counter-rotating streams of particles on stable time-like equatorial circular geodesies 
and that (ii) the inner disc rim is fixed at the smallest possible radius. Keeping the rim at 
the marginally stable orbit of the complete spacetime turned out not to be the correct recipe 
for maintaining the whole disc stable while changing the parameters; namely an instability first 
appears inside the disc, not at the rim. With increasing disc mass, the rim may first shift 
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towards the horizon^] the minimal possible Schwarzschild radius r of the rim goes from 6.00M 
to 3.60M with A4 growing from to 0.42M; the minimal possible circumferential radius R of 
the rim goes from 6.00M to 3.86M with M growing from to 0.29M; and the minimal possible 
proper distance from the horizon d p of the rim goes from 7.19M to 4.26M with M growing from 
to 0.33M. Along the generated sequence of stable discs with minimal inner radii, the radius 
of the marginally bound orbit below the disc was also found to decrease, whereas the photon 
orbit went slowly up; the horizon inflated gradually towards the external source. For A4 greater 
than the values given above, the inner radii of the sequence increase again. When the disc 
mass reaches 1.92M, the position of the inner rim starts to be constrained by the Lagrangian 
point of zero field (below which free circular motion is impossible within a certain radial range); 
this, however, goes up steeply to infinity at A4 = 2M, i.e. the solutions of the given type with 
M > 2M cannot be given the freely counter-rotating interpretation. The Keplerian orbital 
speed was checked to be subluminal everywhere within the resulting sequence of discs, reaching 
roughly 0.7 at its maximum. 

Only the beginning of the above sequence of stable physical discs can be realistic astrophys- 
ically because the mass M. of the actual accretion discs is not supposed to exceed (O.Ol-^O.l)M. 
For A4 between and 0.0722M, the inner rim of our disc sequence follows the marginally stable 
circular geodesic of the corresponding complete spacetime (for Ai > 0.0722M, the rim has to lie 
above the marginally stable orbit to avoid instability of the disc at larger radii). Even with such 
a minor increase of mass, the rim can go down considerably: from r = 6M to r = 4.6162M, 
from R = QM to R = 4.6644M and from d p = 7.1914M to d p = 5.4806M. 

4.2 Oscillations and stability of self-gravitating discs 

The basic criterion for accretion-disc stability, £„ > 0, only concerns perturbations within the 
disc plane. Actually, the frequency k (measured with respect to infinity) of a small free horizontal 
oscillation about a circular equatorial geodesic in a stationary axisymmetric spacetime is given 
by0 



„2v-2X Jlv-2\n 

~ {u*)H 9a<t, ' p 9 Ap ~ (n*)V£ 2 



^-(u')W^p (74) 



([275] and references therein). These oscillations (called epicyclic) have been studied in the 
astrophysical literature since the 1980's, recently as the cause of low-frequency disc modes that 
might explain the quasi-periodic variability observed at some sources probably containing a 



black hole with an accretion disc (see section 2A and the review [157]; the history of frequency 



formulae is described in p75| ). k is zero at radial infinity; in the Newtonian case it increases 
monotonically towards the centre, whereas in relativity it has a maximum at a certain distance 
and then goes back to zero at the marginally stable orbit. Hence, in black-hole sources, the disc 



21 Such a statement is of course coordinate-dependent. Therefore, we also used some physically more relevant 
measures rather then the Schwarzschild radius r only — the equatorial circumferential radius R — pe -I/ ' p ' z ~°- ) = 
re -Vdisc( r -°=K/2) anc j proper radial distance from the horizon (computed along the equatorial plane) d p = 
ff e' A ~ I/ ^ p,z=0 - ) d/9 = f^, r — schw dis,c)( —dr. Note that these transformations are non-trivial as they 

JO f J2M ^/l-2M/r J 

involve the functions fdi ac and Adisc- Consequently, a given r corresponds, for different discs, to different R's and 
dp's (and vice versa). However, the graphs using r and using phy sical radii turned out not to differ much (see the 



direct comparison of all the three coordinates in figure 7 of [274Q. 

22 Note that we consider non-gravitating perturbations, i.e. the metric itself is not perturbed. It would be more 
sophisticated, and also much more difficult, to allow the perturbation to generate its own, self-consistent field. 
However, this would only have a second-order effect on oscillation frequencies. 
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oscillations produced near the inner rim should remain trapped in a region below this maximum 
rather than propagate to the outer parts of the disc fl58|1 . 

In |275[ |, perturbations in vertical direction were found to be also important, in fact, to 
be even more dangerous than the horizontal ones for discs with a greater mass. Perturbations 
perpendicular to the disc plane must be treated carefully because they take place in the vicinity 
of the source which is singular (infinitely thin). Although the symmetry would imply g pU)Z = 
in the equatorial plane, there appears a jump of the field in the normal direction as is usual 
in case of a mass (charge) layer. For the inverted first Morgan-Morgan disc, in particular, one 
finds that (at p > b) 

lim z^ disCi2 = db ^-Jp 2 -b 2 . (75) 

Even in "small" (linear) oscillations the particle spends all of the time outside of the equatorial 
plane where the field is given by ( |75| ) rather than vanishing. The equation of geodesic deviation 
then yields 

- 4(r* tt + r%n)(r* ta + i%,n) - 4(r% + r^o)(r^ + r^o) (76) 

for the "perpendicular" frequency. 

In the static case, the frequency formulae simplify to 

4.V-2X 

k 2 = - (u pp + Apv% - 6»l + 3u p /p) , (77) 



n 



2 e 



P V ,P 
Av-2\ 



1 - pv, P 



4u 2 (l-2pu l 



(78) 



Without the external source, f2j_ reduces to the square of the Schwarzschild orbital frequency, 
f4 = M(M + Vp 2 + M 2 )" 3 , whereas k 2 remains different due to the pericentre precession 
effect. 

Requiring linear stability (also) with respect to vertical perturbations, the range of physi- 
cally acceptable Lemos-Letelier superpositions within the (M, &)-plane is narrowed further (fig- 
ure 1 in |322|). The curve of the marginal vertical stability is more restrictive than that of the 



marginal horizontal stability for the discs with A4 > 0.2296M, it even fully excludes the discs 
with M. > 2M/7. As the boundary of the delimited sector of physical discs (figure 2 in |322j ] ) , 
one obtains a sequence of stable discs whose inner rims lie right on, or very close to, the circular 
geodesies marginally stable with respect to both perturbations. We learned, in particular, that 
for the "realistic" part of this sequence, with relative mass below 0.07, the self-gravity makes 
the oscillations of the inner disc parts faster, and that the region of horizontal mode trapping 
gets somewhat smaller. 

4.3 Redshift from observers in the disc 

The redshift from a given point in the source is often mentioned as a measure of the field 
strength. In Weyl fields, the frequency shift between an observer moving on a circular orbit 
with an angular velocity £1 and an observer at rest at infinity reads 



g = JsE- = e'yji - e -^p 2 n 2 . (79) 

/emit 
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In our case, it is natural to compute the shift from two privileg ed obse rvers within and below 
the disc, the static one (O = 0) and the geodesic one (|fi| = ^ the observer co-rotates 



with the matter in terms of which the disc is interpreted). Figure 3 of |322| shows the curves 



obtained (for discs with different mass and inner radius) from the respective forms of equations 



(7£), g = e u and g = e u \J\S0^i by varying the emitter's orbital radius. 



4.4 Singularity at the rim of the first Morgan-Morgan disc 

The class of Morgan-Morgan [229] static axisymmetric solutions describes the fields of a sequence 



of finite thin counter-rotating discs with densities 

/ \ y / 9 \ "T— 1/2 

w ^ (p < h) = ^±£M(,-^ ( m = i, 2 ....). ( 80 ) 

The zeroth (m = 0) member being clearly singular at the rim, the main attention has been 
devoted to the first (m = 1) member. The latter has been discussed or referred to as a prototype 
of a simple and physically meaningful Weyl field. In addition to the papers mentioned in previous 
sections, the first Morgan- Morgan disc has appeared as the limiting case of a more general class 
of stationary disc solutions [14, 114], while its inverted counterpart has been considered as a 



possible seed for stationary black hole-disc superpositions |p25[| (see the next section). 

The first Morgan-Morgan disc has however turned out to have a curvature singularity at 



the rim [271]. It is caused by too steep a decrease of density (80): the gradient 



w(m) = j^-i)M P f_ pl y^ 

>* 2vr6 4 \ b 2 J { ' 

goes to infinity at p — ► b~ if m = 0, 1. The singularity is inherited by the annular discs, obtained 
by the inversion 

b 2 p ^ b 2 z 

p 2 + z 2 ' p 2 + z l 

(see |322(| ). It was thus recommended [ p71| to consider "higher" (m > 1) members of the 
Morgan-Morgan counter-rotating family in superpositions. More accurately, even these discs 
have certain singularities at their rims — those calculated from higher derivatives of the metric. 
This can already be expected from a "Newtonian look" at densities (80): the m-th member of 
the family has infinite m-th derivative of density at the rim. 

4.5 Inverting the m-th Morgan-Morgan counter-rotating disc 

The potential for the whole family can be most easily written in oblate spheroidal coordinates 
x and y(= cos 6), introduced by 

p 2 = b 2 (x 2 + 1)(1 - y 2 ), z = bxy (0 < x < oo, 1 > y > -1) : (83) 

KA m 

= "X E c£ ) iQ2n(ix)P 2n (y), (84) 

n=0 

where 

(m), (4n + l)(2n)l(2m + l)l(m + n)! 
° 2 " " (n\) 2 (m-n)\(2m + 2n+l)\ {U ~ m} > (85) 
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P2n{y) and Q2n(jx) being the Legendre polynomials and Legendre functions of the second kind, 
respectively; note that it is often suitable to express 

2n ^ 

\Q2ni}x) = P 2 n(»arccot x - i TPk-ifa)P 2n - k (ix). (86) 

k=i K 

Thanks to the invariance of the Laplace equation with respect to the Kelvin transforma- 
tion, one can invert finite thin discs with respect to their outer rims, to obtain annular discs 
(which can in principle be superposed with some central body then). In terms of the oblate 
coordinates, the inversion (|82|) is 

\Jx £ + 1 - y z \Jx l + 1 - y z 

where X, Y are assumed to be related to R, Z in the same manner ( p3| ) as x, y are related to 
p, z. The inverted counterpart of the solution u^^(x,y), satisfying the Laplace equation and 
having the correct asymptotics ~ — j£ , reads 

_ 2 2m + 1 (m!) 2 u^(X,Y) 



n(2m + 1)! y/ x 2 + 1 - y 2 



(88) 



The inverted discs of the Morgan-Morgan class can now be superposed (for example) with a 
Schwarzschild black hole or they can be chosen as seeds in some procedure generating stationary 
generalization of such a superposition (as in the two-soliton inverse-scattering solution described 
in the following section). 



5 Stationary thin discs around rotating black holes 

In order to proceed to stationary superpositions, we employed the Belinskii-Zakharov inverse- 
scattering method [p7| |28| - one of the generating techniques which, in principle, provide 
very general classes of solutions for the problem with two commuting symmetries. In [ |325| , a 
real- two-soliton version of the method was applied to a general Weyl metric (^); the metric 
functions of this "seed" will be denoted by a hat, i.e. v and A, leaving v and A (and u) for 
the resulting potentials. In Boyer-Lindquist-type coordinates (t,r,6,(f>), introduced by (Pq ), the 
obtained metric appears as 

ds 2 = - ^ (Ve^dt + Se~ p d<t)) 2 + (ne dt - Te^dcf)) 2 

+ C 2 e 2X ~ 2i> (|dr 2 + Sdfl 2 ) . (89) 

It involves two functions given by quadratures and it depends on (two) potentials describing 
the seed spacetime and on five independent constants; one can also choose two sign values in 
the derivation. If the seed is flat, one ends up with the Kerr-NUT solution. Restricting to 
reflectionally symmetric, asymptotically flat spacetimes containing a black hole and adjusting 
the coordinates in the most natural manner, the only freedom remains in the seed (y, A) and 
in two constants (one puts C 2 = 1, among others) which we denote M and a as in the Kerr 
solution. The functions present in (]89| ) then read 

A = (r-M) 2 -k 2 =r 2 -2Mr + a 2 , (90) 
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£ = {^rV + ysinhnj + (Ky + asinh vf (=VT + KS), (91) 

V = cosh u - (M/k) sinh u, (92) 

TZ = a cosh (93) 

5 = {l + y 2 )ll + 2aysm\iv -2aVe 20 , (94) 
2a 2 

T = (r 2 - a 2 )P + — (r - M) sinhu + 2aKe 2u , (95) 

ft 



where = V M 2 — a 2 and it and u are given by equations 

A "'' 2A [(r-M)v r sm6 + e y}, (96) 



fc sin 6 A + k 2 sin 2 
u fl 2 



ft sin 6* A + k 2 sin 



2 , 



[(r - M)z> e sin - Av r y] , (97) 



with boundary conditions u(y = ±1) = 0, v(y = ±1) = ±2i>(|y| = 1). 

A number of properties of the above metric were analyzed. The (outer) horizon is given 
by (the larger root of) the equation A = 0. Prom equations of section |3.3| , its area works out at 



A = 4vr (r 2 e - 2 -H(y=D + a 2 e 2 -H(y=D) (98) 
and its surface gravitation and angular velocity with respect to infinity at 

A:7vk 47TCI 

k = —— and wh = — j— cosh 2z>n(y = 1). (99) 
The Gaussian curvature of the horizon reduces to 

(rl - a 2 e 4 ^) [1 - 4i> H H (l)l - 2a 2 

Cn(y = 1) = ^ ^ 4 e 2 -(i) (100) 

(r| + a 2 e 4l/ H(i)) z 

at the axis. This can become negative for rapid rotation, namely for 



e 4£> H (l) J 2 

a>2M v - m l ~^f ] . (101) 

The static limit is located where 

VAP = Usm9. (102) 

Writing the metric in the Weyl-Lewis-Papapetrou form (|]), the asymptotics of the gravi- 
tational potential v and of the dragging angular velocity lj, 

M + M n . 2n 2a(M + 2M) n . 4n . , 

v = + 0{r~ 2 ), u = ^ 1 >- + 0{r- 4 ), 103 

tell us that M + M is the total mass and (M + 2M)a is the total angular momentum of the 
solution; M denotes the mass of the seed spacetime for which we assumed the asymptotics 
v = -M/r + 0{r~ 2 ). 
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In the static case (a = 0), the result represents a superposition of a given seed with a 
Schwarzschild black hole, with u playing the role of an interaction term — u = u + fsdiw; 
A = A + Aschw — u. Thus one can take just the "external source" as a seed, the black hole is 
"supplied" by the soliton method itself. Hopefully, this applies to a general, stationary case as 
well. Actually, our main aim has been to learn whether the Belinskii-Zakharov technique could 
yield the field of a rotating black hole surrounded by an axisymmetric disc. 

It seems that for small enough (though by far not negligible) values of a the solution 
(|8g|)-(|97[) does not contain singularities on or above the horizon. However, it is quite exten- 
sive to calculate and analyse curvature invariants, even if using computer algebra. Coordinates 
were found in which metric is regular on the horizon (they have the meaning of the Kerr ingo- 
ing/outgoing coordinates there), but this is not enough to claim the absence of a singularity. 
Unfortunately, we found the horizon to be pinched rather than smooth in the equatorial plane. 
Indeed, the proper azimuthal circumference 27r(g^)n of the horizon has a sharp minimum in 
the equatorial plane if a > (and M > 0). It is to be clarified whether this pathology is physical 
or the Boyer-Lindquist-type coordinates just do not cover the whole manifold. Most probably, 
however, there appears a supporting surface between the hole and the external disc. 



5.1 



Stationary superposition cultivated from the inverted first Morgan-Morgan 
disc 



Some formulas were already specified to the case of a thin equatorial source in [325], in particular, 
those for the calculation of the energy-momentum tensor and of integral quantities characterizing 
the resulting spacetime. Let us mention below several properties of the metric (|89|)-(|9^) obtained 
for the inverted first Morgan- Morgan disc [275], i.e. for u given by ( |69| ) (with the original 
disc mass now denoted by M). We saw in section 4.4 that this disc is singular at the rim. 



However, its other properties are quite satisfactory and the singularity is only caused by a sharp 
increase of density, so it is reasonable to study first this simple case and only later to embark 
on superpositions generated from "less-singular" seeds of inverted higher Morgan-Morgan discs 
(section |4,5|) . 

On the horizon, 



2M 
7r/c 3 |y| 



(k y + b )arctan 



k\y\ 



bk\y\ 



(104) 



This reduces to i>u(y = 0) = — |^ in the equatorial plane. The main horizon properties 
(9£,99|) only depend on the axial value On(y = 1). In |j273fl it is shown that with growing mass 
or decreasing inner radius of the disc the horizon inflates, its surface gravity weakens and its 
rotation slows down. The presence of the external source weakens the dependence of the horizon 
curvature on a/M. The existence of the sharp equatorial contraction, revealed by the behaviour 
of the proper horizon circumference, has been supported by plotting the isometric embedding 



(H) 



The only undetermined quantities in the metric ( |89|) are the functions u and v. They are 
fixed on the axis by the boundary conditions u(y = ±1) = 0, v(y = ±1) = ±2i>(|y| = 1). On the 
horizon, one easily obtains uu(y) = 2z>h(1) — 2i>n(y), vu_(y) = 2z>h(1) • signy = =Fconst. In order 
that the spacetime be reflectionally symmetric, the seed potential v must be an even function of 
y. Then u is also even in y, whereas v is odd. Details of the courses of u and v are determined, 
according to the equations (96) and (97), by v. In our case, v ( |69| ) is continuous and negative. It 



is maximal on the axis, having z> r > and v t g = there [while v, y (y = ±1) < 0]. In the equatorial 
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plane v is minimal; v r is continuous there, being negative/positive below/above the radius 
p = ^3/26 r = M + Vfc 2 + b 2 ) where v assumes its global minimum of — (2/3) 3 / 2 (M/6). 
The latitudinal gradient is only zero below the inner radius of the disc, while it has a finite 
jump across the disc, given by (|75|) : 

4M6 



lim v g = -(r-M) lim v z = =p— — (r - M)VA - b 2 . (105) 

U->0± ' 2^0± ' ttA z 



At infinity, = —M jr + 0(r 2 ) which implies the asymptotics 

u= 5 + <3(r~ 3 ), u = - + 0{r~ 2 ). (106) 



The values of u and t> at general locations must be computed numerically, one can only 
proceed further analytically in the equatorial plane. There, equations ( p6[ ) and (^) yield lati- 
tudinal derivatives 

2& , 2A „ . . 

u > e = ^7^' = ^7^- 107 

r — M r — M 

Hence, u increases in latitudinal direction from zero on the axis to a maximum in the equatorial 
plane; this maximum is smooth/sharp (ug is zero / has a finite jump there) below/above the 
inner disc radius. The global maximum of u lies where u tT = in the equatorial plane (i.e. 
at p = \/3/2 b; v has a global minimum there). The function v is negative/positive in the 
upper/lower hemisphere, having a finite jump across the equatorial plane; this jump is constant 
- equal to — 4i>n(y = 1) — from the horizon up to the inner disc edge, then it decreases with 
increasing radius. v t g is smooth across the equatorial plane, being positive/negative above/below 
p = 73726. 



The radial derivatives of the equations ([}(]) and (p7[ ) reduce to 
2k 2 

= iTF*V> v ,r = TjV, ( 108 ) 

r — M r — M 

in the equatorial plane. Substituting from (f7l"l), and ( |105|) , one comes to expressions which 
can be integrated to explicit formulae for u(z = 0, p < b), u{z = 0, p > b), v(z = O^p < 
b) = ±2vn(y = 1), v(z = ± ,p > b) (see |7|). While u(z = 0) increases from the horizon 
to a maximum at p(= vA) = \f^>/2 b and then decreases to zero at radial infinity, v(z = 0^) 
increases/decreases to zero monotonically from p = b. 

Using the equatorial expressions for u and v, one can, for example, graph the dependence 
of the equatorial radius of the static limit on M, b and a. We have found that the ergosphere 
widens both with a/M and with M. One can also plot the radial course of the light-like, dragging 
and geodesic angular velocities fimax, u and Q±. 

min 

5.2 Extreme limit 

We assumed an under-extreme case (k > 0) in |325| 1 in deriving the general metric (|89|)-(|97|). 
However, it can be checked where the quantities go in the extreme limit of k — > 0, i.e. a — > M. 
One finds that they behave as in the extreme limit of the pure Kerr solution. 

The event horizon approaches rn = M, its area, surface gravity and angular velocity reach 
the values 

A = 8vrM 2 cosh2i> H , K = 0, w H = ^TJ, (109) 
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where z>h = — for the inverted first Morgan-Morgan disc. The Gaussian curvature of the 



37r6 

horizon works out at 



C H 



4vre 2£>H cosh 3 2z> H (1 + \y\fe, i0u + (1 - |2/|) 3 e" 4pH - 6sin 2 



1 + (\y\ cosh2z>H + sinh 2z>h) z 



(110) 



and reduces to 



C ^ = 1 ) = -2MraS = -T' < 1U » 

at the axis and in the equatorial plane. 

In the extreme limit, u is seen to vanish. The equatorial behaviour of u and v is important 
in discussion of the singularity location S = 0. In [|325| 1 it was shown that this singularity cannot 
lie above the horizon (at least) if a/M is sufficiently small. For our seed (^), however, S(r > rn) 
is positive in general even in the extreme limit. Actually, the first term of (^) then goes over 
to [r — M(r — M)u/k] 2 . Using the expression for the equatorial maximum of u{k — ► 0), the root 
r of the above term must satisfy 

r (M b 2 \ M 

M \ M M 2 ) - M" ^ ' 

This can only hold (for some positive r) if M/M > 3b 2 /M 2 . Such a situation can be considered 
in principle, but it is, in fact, excluded astrophysically: M/M is typically claimed to be much 
less than 1, while a realistic value of 3b 2 /M 2 amounts to the order of ten. 

The static- limit surface ro(9) must be found numerically as equation ( |102[) remains rather 
cumbersome even in the extreme limit. Regarding the limiting forms of u and v, one can, 
however, infer the effect of the external disc on the equatorial value of r$. It is given by 

(r — M)(l — Mu/k) = Mcosh-u, thus r — M > M (114) 

in the extreme limit, which means that rg(y = 0) is greater than its Kerr value (2M). 

It is quite difficult to express, in terms of the parameters M, b, M and a, the contributions 
of the black hole and of the disc to the total mass and angular momentum of the spacetime. 
However, the crucial Komar formula for the angular momentum of the black hole Jh, treated in 



equations (154)-(156) of [325], can at least be integrated in the extreme-limit case, to end up 
with just 

J H = M 2 , M H = M, M = M, J = 2MM. (115) 

Another interesting feature of the k — ► limit of our superposition is the vanishing of 
the "external" gravitational field z> r , A /r on the horizon [|272|| . This effect of expulsion of the 
external (stationary axisymmetric) fields from rotating (and/or charged) black holes, analogous 
to the Meissner effect in (super) conductors, was observed in magnetic fields before (see the 
recent survey 
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6 Concluding remarks 



In sections ^| and || the classes of static/stationary axisymmetric superpositions were presented 
for a black hole with an external thin equatorial disc. They give examples of the effects that 
can occur if the self-gravity of the disc is not completely negligible. The sequence of stable 
static discs with minimal possible radii, generated in section |4.1| , does not (necessarily) describe 
the behaviour of a real accretion disc with its mass increasing (or its inner radius changing). 
Actually, static annular discs need not be of the given (inverted first Morgan-Morgan) type 
(and need not remain within this class when parameters are changed dynamically). Also, real 
accretion discs should be stationary rather than static. The results reviewed in the previous 
sections, however, allow for various generalizations. First, instead of the first Morgan-Morgan 
solution, singular at the rim, one can superpose an inversion of any member of the Morgan- 



Morgan family of static counter-rotating^] finite thin discs (88). Second, a different annular disc 
can be employed, not belonging to the Morgan-Morgan family at all (e.g. the inverted isochrone 
disc of Klein [164]) — and possibly even a thick one (here we only mean the exterior, vacuum 
field). The stationary metric (|89|)-(p7|) applies to any vacuum static axisymmetric seed, too. 

We will try to clarify the nature of the pathology which occurs in the equatorial plane of 
the metric (|8^)~(p7|). If it corresponds to a thin layer of mass rather than to a real singularity, 
it might only be a consequence of the fact that the Boyer-Lindquist coordinates cut off a certain 
part of the manifold. It is more probable, however, that this layer represents a supporting 
surface which bears witness that the system's stationarity is only artificial. What remains an 
open question in any case is the structure of the black-hole interior, altered by the presence of 
the external source. We hope to be able to answer it in future. 

On a general level, one could conclude with the interview with Profs. W. B. Bonnor and R. 
Penrose which Jiff Bicak organized in 1968,0 about mathematics, relativity and "brain drain" 
f30[| (it was a pretty hot topic at the time when "brain strain" was to continue in Czechoslovakia, 
as Prof. Kuchaf could confirm). In his answer to the question of in which direction to expect 
progress in relativity, Prof. Bonnor remarked: "It is also possible that quasars will indeed have 
something to do with gravitational collapse." As I noticed in the literature, at those days this was 
not yet a theme for a gentleman to speak about in a decent astronomical society. Although the 
"collapsed-matter" interpretation of the quasi-stellar radio sources [ |268| | and the "collapsed-star" 
interpretation of X-ray binaries [324 1 were both introduced by "ApJ". (The very importance of 
accretion as an energy source was apparently first recognized by 323 1.) 

Nowadays, ApJ offers "Viewing the shadow of the black hole at the Galactic center" 



within the next few years [105]. Nevertheless, will not a (gravitational) "shine of a black hole" 



be delivered sooner? Gravitational radiation may prove medicinal for rheumatism or elsewhere. 
Will travel agencies cite papers like [292, 175 1 when offering flights intended for "gravitational 
basking"? In any case, the abundance of event horizons in the last years' NASA ADS and 
astro-ph libraries supports Prof. Bonnor's tip. Today it is in fact impossible to follow, at least 
in some detail, but a narrow sector of the black hole research. Not within relativity, but mainly 
not within astrophysics. In the bibliography below, we try to recommend those papers which 
should be interesting both for physicists and astrophysicists, mostly avoiding pure observational 
and pure mathematical references. 

23 However artificial the counter -rot ating interpretation may seem, counter-rotating stellar discs in galaxies have 
become an observational matter [ 304 1 . 

24 It was exactly in 1968 that Prof. J. A. Wheeler coined the word "black hole", but he was only interviewed 
by J. B. several years later [BfJ. 
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By far not all the aspects of ("classical" general relativistic) black holes were mentioned 
in the survey part of this paper. The very question of their origin was only touched in con- 
nection with the cosmic censorship, with the formation of heavy discs and with the decay of 
perturbations. Not a single word about thermodynamics and quantum effects, about cosmology 
and dark matter, about primordial black holes, about economical and ecological aspects, about 
Tunguska impact . . . 

Also, all the black holes considered in the present paper are spheroidal. Recently, toroidal 
holes have been investigated as another intriguing possibility, not only by theoreticians but also 



as a (rather non-orthodox) candidate for active galactic nuclei [284, 254]. Toroidal black holes 
can only exist provided that some of the Hawking & Ellis' assumptions are relaxed. Gannon 
fll §| 1 proved the possibility of a smooth toroidal horizon by dropping stationarity; the period 
for which a toroidal black hole can persist is, however, very limited: the hole in the torus must 



close up before a light ray can pass through 146 ]. Toroidal horizons were really observed as 



a temporary stage of a numerical gravitational collapse ( |295| 1 and references therein). The 
case of a toroidal horizon with a (ring) singularity outside (thus violating the predictability 
assumption) was obtained in |251| 1. Other possibility is to relax the asymptotic flatness and 
allow for a non-zero cosmological constant A. A metric involving (long-term stable) higher- 
genus horizons in an asymptotically anti-de Sitter spacetime is quite simple; it was already 
obtained as a result of a gravitational collapse and generalized to charged and rotating cases 



(see [168 1 and references therein). Astronomical observations indicate that A > 0, but it is not 
excluded that the stabilizing attractive effect similar to that of A < could be provided by a host 
galaxy or by a massive accretion disc around the hole (!). (For many more hints, astrophysical 



connotations and further references, see |284| .) 

In connection with the Weber's putative detection of gravitational waves in 1971, R. 



Penrose speculated about a naked singularity in the Galactic centre |247|| , but at the end he 
stated: "This would remove the observational discrepancies, although it would be a radical 
explanation. . .In comparison black holes can now be regarded as being 'conventional'; indeed, 
for this reason they must be preferred a priori in any attempt at explanation." Our essay has 
been a priori 'conventional', preferring black holes, and of spheroidal topology. Well, Penrose 
then continued: "Yet nature does not always prefer conventional explanations, least of all in 
astronomy." 
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